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Abstract 

We consider a system of nonlinear equations which can be re- 
duced to a degenerate parabolic equation. In the case x G R 2 we 
obtained necessary conditions for the existence of a weakly singular 
solution of heat wave type (codim sing supp = 1) and of vortex type 
(codim sing supp = 2). These conditions have the form of a sequence 
of differential equations and allow one to calculate the dynamics of 
the singularity support. In contrast to the methods used tradition- 
ally for degenerate parabolic equations, our approach is not based on 
comparison theorems. 
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1 Introduction 

We consider the system of equations arising in problems of water purification: 

Oc 

— = DAc- DKdw(cV(divu-3(3c)), (1.1) 

An = (3w-l)/3Vc, (1.2) 

which describes diffusion in media with nonlinear viscosity. Here c > is a 
scalar function, u is a vector, x G Vt C R n , t > 0, D, k, (3 > and uj > 1 are 
several constants. 
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We understand all derivatives in (1.1), (1.2) in the weak sense, consider 
the case n = 2, and obtain necessary conditions for the existence of a solution 
of Eqs. (1.1) and (1.2) with a weak singularity. These conditions, which 
allow us to calculate the dynamics of a singularity, are obtained by using the 
asymptotic expansion of the solution with respect to smoothness. 

Prior to studying (1.1), (1.2), we consider the model example 

du 1 d 2 2 M „. 

* = is?" ■ x 6 (L3) 

One can readily see that a necessary and sufficient condition for the existence 
of a solution finite in x is satisfied for (1.3) (e.g., see [1, 24] and the refer- 
ences therein). On the other hand, from the viewpoint of distributions, only 
the elements of D' that admit multiplication can be singular solutions of a 
nonlinear equation. As is well known [18], in the one-dimensional case this 
restriction leads to an algebra with generators 1, 9(x), eS(x) and x±, A > 0. 
Here 1 stands for all smooth functions, eS(0) = 1 and eS(x) = for x ^ 0, 

( , _ 0, x < x 0, x <0 x \x\ x , x < 

{X) ~ 1, x > 0, X+ ~ x\ x > 0, X -~ 0, x > 0, 

and 9 2 = 9, (x±) 2 = x| A , and 9(x)x\ = x\. 

One can readily verify that both the Heaviside function and e5 cannot 
satisfy Eq. (1.3). Nevertheless, if we set u = (x — a(t))±, then this singularity 
is admissible for A = 1. Such solutions of degenerate parabolic equations are 
well known and are called the heat wave. In particular, the simplest exact 
solution of Eq. (1.3) is the one-parameter family (e.g., see [24]): 



u 



{Qty 1 (x + r]t l ^) + (x -r]t 1/3 )_, r] = const > 0, t > 0. (1.4) 



A natural multidimensional analog of the heat wave (x — a(t)) + is the 
function S + : S + = S for S > and S + = for S < 0, where S = S(x, t) is a 
nondegenerate smooth function. The fact that S is nondegenerate means that 
VS'I^q 7^ for any chosen t, and thus the surface T t = {x, S(x, t) = 0} is of 
codimension 1. Solutions of this type for a wide class of degenerate equations 
have been studied thoroughly (e.g., see [1, 7, 10, 24]). In particular, it is 
known that, starting from some time, the degeneration surface T t becomes 
differentiable [1, 5, 13, 22]. 
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Another possible generalization of (x — a(t)) + to the multidimensional 
case is a singularity of the form \x — a(t)\, x,a e R n . Solutions of such 
form, whose singular support is of codimensions > 1, are called vortex type 
solutions, since they are used for the mathematical description of a typhoon 
eye motion [4] . It is known that vortex type solutions exists for shallow water 
equations [4]. However, as far as we know, the existence problem for such 
solutions of parabolic equations has not been discussed. 

In this paper we obtain necessary conditions under which system (1.1), 
(1.2) have singular solutions of both types mentioned above. 

We study system (1.1), (1.2) by using the asymptotic expansion of the 
solution with respect to smoothness. For linear equations, the idea to expand 
the solution with respect to smoothness was proposed and developed by 
R. Courant, D. Ludwig, V. P. Maslov, V. M. Babich and others. A further 
development of this idea applied to quasilinear equations and special solutions 
lying in D' and admitting multiplication was performed by many authors [7, 
8, 23] starting from Maslov's paper [19]. A possibility to study a more general 
class of nonsmooth solutions of quasilinear equations [2, 3, 6, 9, 11, 12, 21] 
appeared after the construction of algebras of generalized functions, which 
include distributions from D' (V. K. Ivanov, G. Colombeau, Yu. V. Egorov, 
and others). 

A method based on asymptotic expansions of solutions with respect to 
smoothness is very fruitful, since in this method the maximum principle and 
the comparison theorems are not supposed to be valid. On the other hand, 
the following significant drawback is typical: this method implies a system of 
model equations (the so-called Hugoniot type conditions), which is an infinite 
nontriangular sequence in the case of nonlinear equations. 

For example, by setting 

u = ai(t)(x - </?) + + a 2 (t)(x - tp)\ + . . . (1.5) 

in a neighborhood of the singularity support x = <p(t) and substituting (1.5) 
into (1.3), we obtain the system 

(p — — ai, di = 4aia2, a,2 = 3(2<22 + 3aid3), . . . (1.6) 

In the special case a,\(to) = v/(3t^ 3 ), 02(^0) = _ l/(6£o), a«(to) = 0, i > 3, 
for t > t > it follows from (1.6) that 

(p = -r]t 1/s , ai = r?/3t 2/3 , a 2 = -1/6*, a* = 0, % > 3. 
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Thus we arrive at the exact solution (1.4) considered near the left-hand front 

x = —Vjt X l Z . 

In the general case system (1.6) is an infinite sequence of coupled equa- 
tions. Of course, infinite sequences cannot be used for practical computa- 
tions, hence we need to close them. If codim sing supp u = 1, then this can 
be performed by passing to a free boundary problem. In this case conditions 
of the free boundary follow from a local description of the solution near the 
singularity support and from the first Hugoniot type condition. For exam- 
ple, it follows from this procedure that in the region Q t = fit x {t > to}, 
fit = {x,(p-(t) < x < (p+(t)} we need to study Eq. (1.3) with the conditions 

= 0, <p± = 7- 



u 



= u 

x=ip±±Q 



(1.7) 



and a natural initial condition. The first condition in (1.7), which follows 
from the structure of solution (1.5), implies that u and the flow uu x are con- 
tinuous on dQf The second condition in (1.7), which relates the boundary 
velocity to the limit value of the derivative of u along the outward normal u, 
is exactly the first Hugoniot type condition in (1.6). In this statement of 
the problem we can ignore all other conditions in (1.6), since they are au- 
tomatically satisfied due to the fact that the free boundary problem has a 
solution. 

From the ideological viewpoint, this method of closing system (1.6) is sim- 
ilar to closing the Hugoniot type conditions for the shock wave by passing to 
the classical statement of the problem describing the shock wave motion [8]. 
However, here we obtain a problem that differs from the classical one-phase 
Stefan problem (see [20] and [14, 25] for related problems) by degeneration 
of the parabolic operator on the free boundary. 

Another possible approach is to truncate the sequence of Hugoniot type 
conditions. By choosing a number N > and setting all "extra" function 
equal to zero, we can calculate the dynamics of a singularity by using the 
first N Hugoniot type conditions. However, in this cases we need to justify 
this procedure. In general, this problem, which is similar to the well-known 
problem of truncating Bogolyubov-Born-Green-Kirkwood-Yvon chains, re- 
mains open. Nevertheless, the first result in justifying the truncation pro- 
cedure has been obtained [8], and numerical results show a good agreement 
between the exact solution and those obtained by using truncated sequences 
[4, 23] for a reasonably large time interval. 
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In the case codim sing supp u > 1 the truncation procedure is apparently 
the only real method for calculating the singularity dynamics. Indeed, in 
this case, instead of the free boundary problem, we must consider a nonlinear 
equation in a domain with a remote set of codimension > 1. The position 
of this set, on which we pose additional conditions, varies in time and is to 
be determined. As far as we know, such problems have not been studied till 
recently. 

Now let us consider system (1.1), (1.2) and write it in a simple form. To 
this end, we apply (in the D'-sense) the operator div to Eq. (1.2). We assume 
that c has a singularity of an above-listed type and obtain the relation 

A(divu-/3(3w- l)c) = 0, 

which holds for all x including the singularity points. Hence 

divu = j3(3u — l)c + w, 

where w is an arbitrary harmonic function. Now Eq. (1.1) acquires the form 

Oc 

— = D{ div ((1 - k(3{3uj - 4)c)Vc) - «(Vc, Vw)}. (1.8) 

C6 

We assume that oo > 4/3, since only in this case Eq. (1.8) is degenerate and, 
correspondingly, can have singular solutions bounded in the C-norm. 
By changing the variables, we rewrite (1.8) in the compact form 

Oc 

— + {Vw, Vc) = A(c-/ic 2 ), fj, = const > 0. (1.9) 

C6 

Our goal is to derive necessary conditions for the existence of a solution 
with a special local structure. Therefore, we abstract our consideration from 
the initial and boundary (on the outer boundary OQ) conditions. However, it 
is clear that in this case the initial data cannot be arbitrary but must satisfy 
requirement that follow from the Hugoniot type conditions obtained below. 

In the present paper we restrict ourselves to a discussion of Eq. (1.9). 
However, the method proposed here for constructing singular solutions can 
be easily generalized to a wide class of degenerate parabolic equations and 
systems. 
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2 Heat wave dynamics 



In D' it readily follows from Eq. (1.9) that for all t the wave front T t = 
{x,S(x,t) = 0} lies on the degeneration surface E = {(x,y),c = 1/2//} 
corresponding to Eq. (1.9). This fact implies that the asymptotic expansion 
of the solution with respect to smoothness in the neighborhood of T t has the 
form 

c = (2//)" 1 + ai S+ + a 2 S 2 + + ... (2.1) 

Here S G C°° is a function such that VS\ T ^ 0, a± < 0, the coefficients 
t)\ T G C°° have the meaning of the normal derivatives of c on r t , 
and dai/dv\ Tt = 0, where v = — WS/\ VS\ | r is the outward normal. Just 
as in the case of the model equation (1.3), the structure of the solution 
guarantees that the flow (1 — 2[ic)dc/dv is continuous on r t . To find S 
and Oj, we substitute (2.1) into (1.9) and expand Vw in the Taylor series on 

r t : 

Vw = V Q + V 1 S + V 2 S 2 + ... 

Now by setting the coefficients of 9(S), S + , S+, . . . equal to zero, we obtain 
the infinite nontriangular system 

{S t + (V , VS) + 2/i|V5|V}| 5=0 = 0, (2.2) 
K + (V , V^)a J? + 2(i + 1) 2 //| VS\ 2 ai a i+1 + F t } = 0, % = 1, 2, . . . 

Here ^ = ip(x,t) is a parameter along T t that is determined in each chart as 
the solution of the equation 5 = 0, and Fi are functions of S, a±, . . . , and 
their derivatives. In particular, we have 

F l = a 1 (V 1 ,VS) + 2/ia 2 1 AS, 

F 2 = 2a 2 (V 1 ,VS) + a^(V 1 ,V^)+a 1 (V 2 ,VS) 

+fi{Q ai a 2 AS + \2a\\VS\ 2 + |V^| 2 (a?)^ + A^(a 2 ) J. 

Equations (2.2) are written in the form that is inconvenient for analysis 
and numerical calculations, since in (2.2) we need to calculate the trace on 
r t . This can be avoided by taking into account the geometry of the front 
r t . For example, if for t G [0,T) the front can be uniquely projected on the 
a; 2 -axis, then we can set S — X\ — ip(x 2 ,t), obtain VS = (1, —(p X2 ) T , £ = x 2 , 
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dj = a,i(x2,t), and rewrite Eq. (2.2) as 



<Pt + w X2 
a lt + w X2 



xi=<p 



ip X2 = 2//(l + (pi )ai + w. 



xi=ip 



(2.3) 



xi=ip 



xi=¥)' 



+2/i( - ai^^ - 2(a\) X2 ip X2 + 4(1 + <p\ )aia 2 ) = 0, . . . , 



If r t is a closed smooth curve without self-intersections, then we can avoid 
a local description of the front by setting S = t — $(x). Then = aj(x) and, 
instead of (2.2), we obtain 



1 - (V$,Vw> 
-(Vai, Vw) 



t=<f> 



+ 2/iai|V$| 2 = 0, 
■ai<Vi,V$) 



(2.4) 



+2/iai{4(V$, Vai) + aiA$ - 4a 2 |V$| 2 } = 0, . . . , 



where V l = (S7w t )\ t=q ,. 

By writing the Hugoniot conditions in the form (2.3) or (2.4), we arrive 
at more trivial problem of posing the initial conditions. Namely, for (2.3) we 
specify the values of ip±, ai, . . . for t — 0, and for (2.4) we specify the values 
$ = and ai — a®, . . . on the initial curve T whose position is assumed to 
be known in advance. 

Thus we arrive at the statement. 

Theorem 1 For the existence of a heat wave type solution of Eq. (1.9), is 
necessary that the Hugoniot type conditions (2.2) be satisfied on the front T t . 

Note that ai|VS'|| rt = —dc/du\ T and v v = S' t /|V5'| is the velocity of the 
boundary r t along its outward normal. Therefore, the first Hugoniot type 
condition can be written in the form 



dc 



dw 
dv 



(2.5) 



As was already noted in Introduction, condition (2.5), together with the 
condition c\ v = (2fi)~ 1 , allows us to reduce calculations of the heat wave to a 
one-phase free boundary problem. Obviously, the solvability of this problem 
implies that all Hugoniot type conditions are satisfied. 
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Remark If instead of (1.9) we consider the equation 



! + <v,vo 



(V,*(c)Vc>, 



where k = A<)C 7 (1 + 0(c)) as c — > 0, k > 0, and 7 > 0, then by analogy with 
(2.1) the asymptotic expansion of the solution with respect to smoothness 
has the form 

C = CLiSl 



+ a 2 ST 1+1 + ... 



In this case the first Hugoniot type condition is similar to (2.2), namely, 



{ St + (V,VS)-^\VS\ 2 aj} 



5=0 



0. 



However, now the closure of the sequence of Hugoniot type conditions 
leads to a one-phase Stefan type problem, which consists of the original 
equation considered in the domain Q t = Q t x {t > t }, Q t = {x, S(x,t) > 0} 
with the initial condition and the boundary conditions on the free boundary 
r t = dtt t : 

0, v„ = (V,v) t ' - 1 ' 



ov 



where /| r is understood, just as in (2.2), as the limit obtained by passing to 
T t along the inward normal. Obviously, we can rewrite the last term in the 
Stefan condition in the form — k[c)c~ x dcj dv\ Vt . 



3 Dynamics of a vortex type singularity 



The solution of Eq. (1.9) with a vortex type singularity can be written in the 
form 

C = C°(x,t) + yj 'S(x,t)c l (x,t), (3.1) 

where S, c°, c 1 are some smooth function. We assume that for each fixed t 



S > 0, VS 



5=0 



0, 



S" 



5=0 



> 0, 



(3.2) 



and moreover, T = \J t T t C R+ x R 2 is a curve in the C°°-class. Here T t = 
{x, S(x, t) = 0} is the singular support of the solution and S" is the Hessian 
of S. By a(t) we denote a vector-function, which describes the position of 
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a singularity at each time instant t, i.e., we have T t = {x,x + a(t) = 0}. 
We perform the change x' = x + a(t). Then, omitting the prime on the new 
variable, instead of Eq. (1.9), we obtain the equation 

c + (d + Vw, Vc) = A(c-/ic 2 ). (3.3) 

Here and in the following, the dot indicates the derivative with respect to t. 

Substituting (3.1) into (3.3) and grouping smooth and singular terms, we 
obtain the relation 

D° + S'^D 1 = 0, (3.4) 

where 

D° = c° + (V, Vc°) - AG , V = a + Vw, 

D l = -kG 1 +AS(^(VS, V)c l - (V^VG 1 )) (3.5) 

+ AS 2 (-AG l + c 1 + (V, Vc 1 )), G° = c°- /i(c° 2 + Sc l2 ), 
G 1 = (l-2fic°y, k = 2SAS - \VS\ 2 . 

The smooth functions D l satisfy relation (3.4) at the point x = only if for 
any N > we have 

Qk 



D° 

dx a 



x=0 



0, \a\ — k, k < N, (3.6) 



dx a 



= 0, |a| =k, k < N, (3.7) 

x=0 



where, as usual, a = (0:1,0:2) is a multiindex. 

We write the functions S and c\ % = 0, 1, as the formal Taylor expansions 

00 00 

e = xy*(M)> s = E%*). (3-8) 

fc=0 fc=2 

pi = S c «(^ Q . ^ = E 

|a|=fe |a|=fe 

In view of conditions (3.2), we assume that sn = 0, S20 > 0, and S02 > 0. 

In what follows, we essentially use the fact that the vector-function V is 
of a special structure. Namely, since w is a harmonic function, we have 

00 

w = Y,W k (x,t), W k = u a (t)x a , (3.9) 

k=0 \a\=k 
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where Wk are harmonic polynomials. Therefore, in particular, Wq = Wo(t) 
and Wi(x,t) are arbitrary functions, and 

W 2 = uj 2 o(xI - x\) + U11X1X2, (3.10) 
W 3 = uj^ xi(x\ - 3x1) + ^03^2(^2 - 3x1), 

where u) a = oj a {t) are arbitrary smooth functions. Thus in the Taylor expan- 
sion 

00 

V = Y,V k (x,t), V k = Yl v°(t)x a , (3.11) 

fc=0 \a\=k 

we must take into account that Vq = a + VWi and Vi = VWj+i, % > 1. 

Substituting expansions (3.8)-(3.11) into (3.6) and (3.7) and setting the 
coefficients of x a equal to zero, we obtain an infinite system of equations. One 
can readily see that this system is not triangular, i.e., the first N equations 
contain M > N unknowns. The inequality M > N is necessary for the 
existence (at least, for small t) of a solution of these equations, but, in general, 
is not sufficient. The matter is that the compatibility conditions for these 
equations contain less than M unknown functions. Since we do not know in 
advance whether the equations obtained are solvable, a detailed analysis of 
these equations is required. 

To avoid too cumbersome formulas, we shall successively analyze the 
equations obtained for k — 0, 1, . . . 

We consider (3.7) for k — 2. Since 

00 

k — ^ Kk(x, t), k 2 = 2S 2 AS 2 - \VS 2 \ 2 = 4s 2 ()So2^ 2 , (3-12) 

fc=2 

we have the relation 

S20S02Cq(1 - 2/ICq) = 0. 

In view of (3.2) and the natural assumptions that 

cl ± (3.13) 

(otherwise, the singularity of the solution (3.1) becomes weaker that \x\), we 
derive the relation 

eg = 1/2//, (3.14) 
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which readily implies that the condition /i > is necessary for the existence 
of a physically meaningful (c > 0) solution of the form (3.1). 

Now we consider (3.7) for k — 3. Since Gq = in view of (3.14), we have 
the relation 

k 2 G{ + AS 2 (VS 2 , X7G\ - l -V cl) = 0. (3.15) 

Here and in the following, relations of the form Pk(x, t) = J2\ a \=k Pa(t)x a = 
mean that the coefficients P a (t), \a\ = k, are equal to zero. We also use 
expansions of the functions G l that are similar to (3.8): 

G l = EQ(x,t), G\ = J2 9a(t)x a , i = l,2. (3.16) 

fc=0 \a\=k 

In particular, by taking into account (3.14), we have 

Gl = -^clpl+pf+pf), Gl = 0, (3.17) 
G\ = Gh-2M« + p|p!). 

By substituting explicit formulas for k 2 , S 2 , G\ into (3.15) and matching the 
coefficients of like powers of x, we obtain the four equations 

6/ic° + V , a = 0, 2/i( So2 / S20 )^c° + Vb.a + 4/ic° = 0, (3.18) 

where a attains the values (1,0) and (0,1), V 0ta = Vq,i an d 9a — 1 f° r 

a = (1,0), Vo, a = Vo,2 and g a = —1 for a = (0, 1), and V 0ti are components 
of the vector V = (V 0; i, V 0t2 ). 

Obviously, Eqs. (3.18) are compatible only for s 2 o = so2- Without loss of 
generality, we obtain 

S20 = s 02 = I. (3.19) 

The (3.18) implies the relation 

V + 6^Vp? = (3.20) 

between the velocity a of the singularity and the coefficients in expansions 
(3.8) and (3.9). 

Next, from (3.6) with k = we obtain 

eg + <V , Vrf) = AG° 
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In view of (3.14), (3.19), and (3.20), this immediately implies 

c l 2 = |Vp?| 2 . (3.21) 

Now let us consider (3.7) for k — 4. In view of (3.19), we have 

x 2 R 2 + S 3 Li = (3.22) 

with homogeneous polynomials R2 and L\ of degree 2 and 1, respectively. 
One can readily see that the compatibility condition for the equations ob- 
tained from (3.22) is 

AZi = 0, 

where A is the 2x2 matrix with coefficients A n = — A 2 2 = S21 — s 03) 
A12 = A 2 i = S30 — S12 and \ = (ho,loi), where l a are coefficients in the 
polynomial L\. The choice l\ — implies that the solution constructed is 
trivial. Therefore, det A = 0, and hence, in turn, 

S 3 = x 2 Q 1 (x,t), (3.23) 

where Q\ is a homogeneous polynomial of the first order. 
By using (3.23) and the obvious identity 

(x, Vp n ) = n Pn , (3.24) 

we transform (3.22) as follows: 

5p°2 + Q1P1 + + ~W 2 = \x 2 ((VQ^ Vp?) + -±j9 ), (3.25) 

where 

d = -c 1 + Q f i(\7p° 1 ,Vp 1 1 ) + AG 1 2 . 

Relation (3.25) implies three equations for the coefficients in expansions (3.8). 
We combine these equations and, in particular, obtain the system 

AVQ! = -^AVp\-f, (3.26) 

where A is the 2x2 matrix with coefficients An = A22 = Cq 1; ^.12 = — A21 = 
c io> / = (/i> ~h), h = 5c? x + wu/n, and f 2 = 5(c° 20 - cg 2 ) + 2w 20 //i- 
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By using (3.21) and the assumption (3.13), we find from (3.26) the coef- 
ficients in the polynomial Q\\ 

VQi = -\vp\-^ ^o = A~ 1 f. (3.27) 

Now the third equation obtained from (3.25) can be reduced to the form 

cl + ^/icjApQ = 0. (3.28) 

Next, we consider (3.6) for k — 1 and arrive at the relation 

p° - 2/i(Vrf, Vp° 2 ) + (Vi, Vp?) + 2MAp° 2 - 4^) = 0, (3.29) 

where is a homogeneous polynomial such that V^i = c^ipo. One can 
readily see that the system of relations (3.21), (3.28), (3.29) is consistent 
only if the compatibility condition 

(Vp?,£Vp?>=0 (3.30) 

is satisfied. Condition (3.30) can be readily obtained by differentiating (3.21) 
with respect to t and calculating Vp? with the help of (3.29). In (3.30) 

£ = 2^ + ^ + 87-^', 

/ is the unit 2x2 matrix, the double prime indicates the Hessian of a 
function, and / is the 2x2 matrix with coefficients f u = — f 22 = /2 and 

/12 = /21 = fl- 
it follows from (3.30) that 

BVp\ = 7(jVp? ± , (3.31) 

where f L = Tf is the vector orthogonal to /, T is the rotation matrix, i.e., 
T11 = T 2 2 = and T i2 = — T 2 i = —1, and a = a(t) is an arbitrary scalar. 

Let us transform (3.31). To this end, we note that B can be rewritten 
in the form B = 7(B° + /i -1 ^'); where B° is the matrix with coefficients 

*->ll ~~ ,c 20 dc 02) ^22 — UL 20 ^ ' c 02> °12 — u 21 ~ uc ll- 

Moreover, the following relations hold: W 2 Vpi = DuJ 2 and i3 Vp? = 
B 1 C 2 + c 02 r, where w 2 = (u 20 ,un) T , C 2 = (c%,c° u ) T , and r = (-5c? , 7c% 1 ) T , 
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D and B 1 are the matrices with coefficients Dn = 2D 22 = 2c5 , -D21 — 
-2D 12 = -2c° 01 , B\J7 = B] 2 /Q = c? , and Bj 2 /6 = -B\Jh = c° 01 . Therefore, 
we can rewrite (3.31) in the form 

DoJ 2 = n(aT - B Q )Vp\ = niaVpl 1 - - B 1 C 2 - c° 02 r}. (3.32) 

Both matrices D and B 1 are nondegenerate in view of the assumption (3.13) 
(det-D = 2cJ 2 and det-B 1 = 6(5cq 2 + 2c5 2 ) = f Qd). However, we must take 
into account that any condition imposed on W 2 implies restrictions on the 
boundary values of u. So, avoiding the appearance of such restrictions and 
inverting B 1 , we obtain from (3.30) the algebraic relations 

c° = i{-2ac? c[li + c° 2 (5c; o 2 + 7c° 1 2 )-2c^ 2 o//i}, (3.33) 

c° n = -L{- a (5c° 01 2 - 7c% 2 ) - 2Ac 02 c° 10 c° m + Ac° 10 c° 01 u; 20 / pi 

-(7c° 10 2 + 5c° 01 2 )ujn/ri. 

Obviously, since the higher-order terms of the expansion contain differential 
equations for the coefficients of p®, we must prove the validity of relations 
(3.33). 

First, we note that (3.32) allows us to simplify Eq. (3.29). Namely, we 
invert the matrix D, calculate V\ and ipi, and, according to (3.21), set 

c io = Cosffiy 9 ) c oi — cl cos (f, (3.34) 

where tp — tp(t) is a new unknown function. Then after some calculations we 
transform (3.29) to the final form 

ip + 7//0- = 0. (3.35) 

To prove that (3.33) is consistent, we consider (3.6) for k = 2, which leads 
to the equation 

P°2 + (Vo, Vp°) + (V ± , Vp° 2 ) + {V 2 , Vp?) = AG", (3.36) 

and (3.7) for k — 5, which leads to the relation 

6fju%S4p1 + x 2 TZ 3 = 0, (3.37) 
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where 1Z 3 is a homogeneous polynomial of the third order. By analogy with 
(3.23), we can readily verify that (3.37) implies the relation 

S 4 = x 2 Q 2 (x,t), (3.38) 

where Q 2 is a homogeneous polynomial of the second order. Then (3.37) 
acquires the form 

- 4cSp?Q 2 + TZ 3 + x 2 {cl((VQ 2 , Vp°) - p?AQ 2 ) - iT^X + A} = 0, (3.39) 
where 

n 3 = -2{4py 2 + 7cy 3 + 7 P ipi + QM^p\ + clQi) + iicy 2 ) 

+(W 2 {p\ + 2Q 1 c 1 Q ) + 3W 3 c 1 /2)/^}, 
£1 = Qi(cl(VQ 1 ,V P ° 1 )+ 29 / f i) + \VQ 1 \ 2 cy i /2 
+ ^(VQ 1 ,VGl-(p\V + c 1 V 1 )/2) 
+ A r 1 (AG*- (V ,Vpl) - (VuVpD). 

Relation (3.39) is equivalent to four scalar equations. The compatibility 
condition for these equations is 

Aq 2 = -g 02 Vrf ± - F/Acl 

where A is the same matrix as in (3.26), q 2 = (Q2o,Qu) T , Q a , \ a \ — 2, are 
the coefficients in the polynomial Q 2 , and F is the vector with components 
F\ = T03 — r 2i an d F 2 = r 30 — r 12 , where r a are coefficients in 1Z 3 . Hence we 
find 

Q20 = Q02 + (F , Vpf)/^ 3 , q u = —{F, Vp?)/4^ 3 , (3.40) 

and see that qo 2 is arbitrary. 

Now it follows from (3.40) that the other two relations that follow from 
(3.39) are 

Vp\ =nF- 6/^02 Vp°. (3.41) 
Here T is the vector with components 

•Fi = r 30 + l l 10 - (FAc%Vp\ L + c° 01 Vp° 1 )/Ac 1 \ 
T 2 = r 3 + ^ 1 -(F,2c° 1 Vrf ± + c? Vrf)/4ci 2 , 
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where l\ are coefficients in C\. 

In turn, relations (3.40) allow us to rewrite Eq. (3.36) in the form 

Pi + M + i6g 02 4 V = 0, (3.42) 

where the coefficients in the polynomial M are independent of go2- 

We differentiate relation (3.32) with respect to t and calculate the deriva- 
tives pi and p 2 according to (3.28), (3.34), (3.35), and (3.42). After elemen- 
tary calculations, we obtain 

(& - 7^a/\pl)Vp^ + (32cJ 2 go2 - 14/^ 2 )Vp? = T 1 . (3.43) 

Here T x = [i~ x Dlo 2 — AiVpi and Ai is the matrix with coefficients Ain = 
7M 20 - 5M 02 , M 22 = 7M 02 - 5M 20 , and Mi 2 = M 2 i = 6M 1U where M a are 
coefficients in the polynomial M. 

Obviously, (3.43) implies the differential equation 

&-7v U TA& = (F 1 ,Vp* L )/<g (3.44) 

and the relation 

q 02 = ( 7 ^ 2 + (F\ Vp\)/2cf)/lQcf, (3.45) 

which guarantee that (3.33) and (3.36) are compatible. 

In turn, relations (3.33) and (3.45) mean that (3.42) can readily be re- 
duced to the scalar equation 

4 + 7fia 2 + M 02 + (F\ Vp\)/2cf = 0. (3.46) 

Let us summarize the preceding. Relations (3.14) and (3.20) mean that 
Eq. (1.9) degenerates on the singularity support (since we have c = c° = 1/2// 
there). The construction of the first terms Pi,p 2 and c} j1 p\ in the asymptotic 
expansion of the solution with respect to smoothness (3.1), (3.8) is reduced 
to solving ordinary differential equations (3.28), (3.35), (3.41), (3.44), and 
(3.46) with regard to algebraic relations (3.21), (3.27), (3.33), (3.40), and 
(3.45). 

We introduce the notation y — (y 1 , . . . , y 4 ), where 

yi = a(t), y 2 = c\ {t), y 3 = c 1 01 {t), y 4 = c ° 02 (t), 
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and z = (z\, . . . , z-j), where Zj is a function from c^,(t), \a\ = 3, or from c^(t), 
\f3 1 =2. Then we can rewrite the above equations as follows: 

di = — 6/zCq siny? — u w , 
d 2 = — 6/iCq cos ip — uqi, 

= -7/iyi, (3.47) 
9 

Co = iic 1 —(y 1 sm2(p-12y 4 + 2u 2 o/n), 

Vi = Vi (y, ^20, ;</?,/*) + 4^(^30,^03;^/^) 

+ Mj(cj 2 o, ^ii; V 5 , A*) + CoiVi(z; </?, //), i = 1, . . . , 4. 

Here d' = cos2y? — 6, Pj(r; </?,//) are homogeneous second-order polynomi- 
als in r, Li(£;ip, [/,), Mj(£; </?, //), and Ni(z;(p, (/,) are homogeneous first-order 
polynomials in £ and z. The coefficients of these polynomials are uniformuly 
bounded smooth functions of </? and the parameter \i > const > 0. 

The right-hand sides in (3.47) can readily be calculated by using comput- 
ers and "Mathematica" software. Nevertheless, the explicit formulas thus 
obtained are too cumbersome, and we do not write them here. We also note 
that system (3.47) is nonclosed, since it contains coefficients of the polyno- 
mials and p\. 

Finally, it was proved that the function S is of the form S = x 2 (l + Q\ + 
Qi + . . .), and thus implies that the assumptions (3.2) are satisfied. 

We consider the higher-order terms in the expansion with respect to 
smoothness following the above scheme and arrive at the statement of the 
theorem. 

Theorem 2 For the existence of a vortex type solution of Eq. (1.9), it is 
necessary that the singularity support x = —a(t) be the curve of degeneration 
of Eq. (1.9), i.e., c\ x= _ c ^ = 1/2//, and that the limit values of the derivatives 
of c as x — > —a(t) satisfy Hugoniot type conditions the first of which have 
the form (3.47). 

For the initial data, it should be noted that relations (3.21) and (3.33) 
impose restrictions on the initial values of c° , \a\ — 1, 2, while | t=0 , \a\ < 1, 
and uj a (t), t > 0, can be chosen arbitrarily. 

The plots show the trajectories (corresponding to different initial data) 
of the singularity support motion calculated without regard for the drift (for 
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ui a — 0). These plots were calculated by using the truncated system (3.47), 
where all "extra" functions c°, |a| = 3, and cL \j3\ = 2, were set to be equal 
to zero. 

The author is grateful to A. I. Koshelev and V. A. Galaktionov who 
stimulated his interest in this problem. 
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